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Abstract 

The radio frequency 1.4 GHz transition of the atomic hydrogen is one of the impor¬ 
tant tracers of the diffuse neutral interstellar medium. Radio astronomical observations 
of this transition, using either a single dish telescope or an array interferometer, reveal 
different properties of the interstellar medium. Such observations are particularly use¬ 
ful to study the multiphase nature and turbulence in the interstellar gas. Observations 
with multiple radio telescopes have recently been used to study these two closely re¬ 
lated aspects in greater detail. Using various observational techniques, the density 
and the velocity fluctuations in the Galactic interstellar medium was found to have 
a Kolmogorov-like power law power spectra. The observed power law scaling of the 
turbulent velocity dispersion with the length scale can be used to derive the true temper¬ 
ature distribution of the medium. Observations from a large ongoing atomic hydrogen 
absorption line survey have also been used to study the distribution of gas at different 
temperature. The thermal steady state model predicts that the multiphase neutral gas 
will exist in cold and warm phase with temperature below 200 K and above 5000 K 
respectively. However, these observations clearly show the presence of a large fraction 
of gas in the intermediate unstable phase. These results raise serious doubt about the 
validity of the standard model, and highlight the necessity of alternative theoretical 
models. Interestingly, numerical simulations suggest that some of the observational 
results can be explained consistently by including the effects of turbulence in the mod¬ 
els of the multiphase medium. This review article presents a brief outline of some of 
the basic ideas of radio astronomical observations and data analysis, summarizes the 
results from recent observations, and discusses possible implications of the results. 

Key words: ISM: atoms - ISM: general - ISM: structures - radio lines: ISM - turbulence 


1 


Introduction 


The interstellar space in a galaxy is far from empty. Even if the density is very low 
(~ 10 -4 — 10 6 cm -3 ), the region between the stars is filled with the interstellar medium 
(ISM) consisting of gas, dust, charged particles and magnetic fields. Astronomers inferred 
the existence of interstellar gas after the discovery of stationary absorption lines of ionized 
calcium (Ca II) towards the spectroscopic binary <5 Orionis (IHartmanni 19041) . Any spectral 
line from stars of a spectroscopic binary system should have a periodic wavelength shift due 
to the Doppler effect for the motion of the stars. The stationary Ca II absorption line was 
hence taken as an indication of interstellar gas. Edward Emerson Barnard first carried out 
a systematic study of the Galactic ISM, and proposed the existence of “intervening opaque 
masses” between the stars to explain “dark nebula”. He used his expertise of astrophotog- 
raphy to produce the first images of dark nebulae and published the first catalogue of such 
dark clouds iBarnardl . Il919 ). Over time, further imaging and spectroscopic observations, 
made the existence of widespread ISM, including gas and dust clouds, quite evident. 

It is now well established that the process of star formation through the collapse of pro- 
tostellar clouds is never totally efficient. This leads to the existence of residual gas around 
the stars. The radiation and mechanical energy inputs from the stars, in turn, influence 
the properties of the ISM. Processes like stellar winds, episodic ejection of mass during the 
stellar evolution and supernova explosions transfer material and energy from the stars to the 
ISM. Similarly, ultraviolet radiation from hot, young stars and cosmic rays heat and ionize 
the ISM. Due to all these processes, the ISM is strongly coupled to the stars in a galaxy, 
and play a crucial role between the stellar and the galactic scales. Hence, understanding the 
properties of the ISM is very important from many considerations in astrophysics. 

Multiphase nature and turbulence are two of the key ingredients of the ISM physics. In 
the standard model, multiple phases of the ISM, with different densities, temperatures and 
ionization states, coexist in rough thermal pressure equilibrium ( Field . 19651: Field et al. . 


19 6 9). In a multiphase medium, the cold dense neutral gas (cold neutral medium; CNM) 


is embedded in either neutral or ionized low density warm gas (warm neutral and ionized 
medium; WNM/WIM), which again may be embedded in much lower density hot ion¬ 
ized medium (HIM). Considering different cooling and heating processes, it can be shown 
that the neutral atomic hydrogen gas (H I) can be in thermal steady state in one of the 
two stable ranges of temperature: (i) ~ 40 — 200 K for CNM and (ii) ~ 5000 — 8000 K 
for WNM. The H I at intermediate temperatures is unstable, and, due to runaway heating 
or cooling, moves to either cold or warm phase dMcKee & Qstriken. 1 1977c IWolfire et all 
1995L120031) . The ISM (of the Milky Way as well as of other galaxies) is also known to 


have a c l umpy , self-similar, hierarchical structure over several orders of magnitude in scale 
(ILarsonL Il98ll: iFalgarone et all 1 19921: iHeithausen et all 119981) . Various d irect and indirect 


observations sugg est that these structure extends down to < 10 AU scales (IHill et. a1.L 12005 


Dutta et all 12014). On the theoretical front, the ISM is known to be turbulent, and hence 


is expect ed to exhibit such density s tructures and velocity fluctuations oyer a wide range 
of scales ( Crovisier & Dickev . 1983 : Greenl 1993 : Lazarian & Pogosvan . 2000l) . However, 
the nature of the turbulence in different phases of the ISM is not yet understood properly. 

This article presents an overview of some recent results, from radio astronomical obser¬ 
vations, which have important implications for our current understanding of the multiphase 
turbulent ISM. A significant fraction of the observations, on which these results are based, 
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Figure 1: A toy model of atmospheric opacity to electromagnetic radiation at different wavelength. 
Note the low opacity for optical (A ~ 400 — 700 nm) and radio waves (few cm to ~ 10 m). 


were carried out with the Giant Metrewave Radio Telescope (GMRT; LSwarup ct ah, 199 lh . 
a world-class radio telescope located near Pune, Maharashtra, in India. A few important 
aspects of the observational techniques in radio astronomy, and a brief overview of the 
GMRT are outlined in the “Observation and Analysis” section. The results, related to the 
ISM turbulence and the multiphase nature, are summarized in sections titled “Turbulence 
in the diffuse neutral ISM” and “Temperature of the diffuse neutral ISM” respectively. Fi¬ 
nally, a discussion on the broader aspects of these results, and on possible future directions 
are presented in the “Discussions and Conclusions” section. 

Observation and Analysis 

Historically, astronomy was confined to observations of the sky mostly with optical tele¬ 
scopes, probing only a narrow window of visible lights of the whole electromagnetic spec¬ 
trum. With technological progress and developments, today astronomers are, however, ca¬ 
pable of observing, using ground- and space-based telescopes, across the whole spectrum 
- 7 -rays, X-rays, ultraviolet, visible, infrared, microwave and radio emission originating 
in different astrophysical processes from a variety of sources. Due to low or negligible 
atmospheric opacity at the radio “window” (A ~ a few cm to about 10 m; see Figure |T|), 
observations at these wavelengths are possible to carry out from ground-based telescopes. 

Radio astronomy emerged as a subject in the first half of the last century, with the pio¬ 
neering work of Karl Guthe Jansky, who first reported the detection of celestial radio signal 
(Jansky. 1933). His contribution has been recognised in many ways - from the unit of flux 
density Jansky (1 Jy = 10" 26 W m 2 Hz ' ), to recent renaming of the upgraded radio 
telescope as the Karl G. Jansky Very Large Array (VLA). Over time, radio observations 
have revealed thermal and non-thermal emission from known as well as new classes of 
sources including the Sun, supernova remnants, pulsars, regions of ionized hydrogen (H II 
regions) associated with star formation, active galactic nuclei, radio galaxies and clusters 
of galaxies. To_get an overview of the exciting results from radio astronomical studies, see 
Burke & Graham-Smithl (2014). For the purpose of this article, let us now consider a very 
important radio frequency transition, the famous H I 21 cm line, in greater details. 


Atomic hydrogen 21 cm radiation 

In the radio frequencies, the most useful tracer of the ISM is the H I 1.4 GHz (21 cm) 
line. This is not only because hydrogen is the most abundant element in the ISM, but also 
because it is possible to extract a lot of useful information (e.g. temperature, column den- 
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sity, velocity, magnetic field strength) about the physical properties of the ISM from H I 
21 cm observations. In 1944, Hendrik van de Hulst predicted the possibility of detecting 
H I 21 cm radia t ion of celestial origin (Ivan de Hulstl 1 1945b . and it was first observed by 
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Ewen & Purcelll (1195 lb . This line emission (or absorption) is caused by the transition be¬ 
tween the two hyperfine states of the l 2 ,5'i ground state of hydrogen. These two states with 
parallel and anti-parallel electron/proton spin configuration (i.e. total spin angular momen¬ 
tum F — 1 and 0) have slightly different energy (see Figure [2]). The energy difference hv 
corresponds to the frequency of the transition v w = 1420.405752 MHz (~ 21.1 cm). 

The 21 cm transition of neutral hydro¬ 
gen has an extremely small probability of 
2.!) x 10 1 \s 1 . So, the natural lifetime of 
the excited state for a single isolated atom 
is about 10 million years. But, because 
of the processes like collisions and inter¬ 
action with the background radiation field, 
the lifetime can be considerably shortened. 

This, and the fact that the total number of 
atoms in the ISM is quite large, make the 
21 cm emission (or absorption against a 
background source) readily detectable. 

In thermodynamic equilibrium, the 
number density of atoms in the up¬ 
per and lower energy state are related 
by the Boltzmann distribution, n i = 

37i 0 exp(— hvi 0 /kT). Even if the gas is not 
in thermodynamic equilibrium, one can 

define a characteristic temperature, called the spin temperature (T s ), that satisfies the ob¬ 
served distribution between the two spin states. 

Either collision or radiative mechanisms may cause direct transition between the two 
hyperfine states. Alternatively, higher energy levels of hydrogen atom can be excited from 
one of the hyperfine states, followed by de-excitation to the other hyperfine state via Ly-o: 
or higher Lyman lines transitions (see Figure [2]). The probability of transition via Ly-a pho¬ 
tons depends on the intensity profile of radiation near the Ly-a frequency. Radiation which 
has been scattered many times in the ISM will have a profile dependent on the kinetic tem¬ 
perature (Tfc) of the gas. This may tend to couple T s to T k , particularly in low density gas. 
In high density regions, the collisional mechanism dominates, and may similarly couple 
T s to Tfc dFieldl. Il958l) . W hether or not T, is same a s T k , for a constant T s along the line 
of sight, it can be shown ( Kulkarni & Heilesl 1988 ) that T B (v) = T s [ 1 — e~ T ^], where 
Tfj(iy) = l{v)c 2 /2kv 2 is the observed brightness temperature, r{v) is the optical depth at 
frequency u, I(v) is the intensity, k is the Boltzmann constant and c is the speed of light. 
At low optical depth limit, the spin temperature will be T s (u) ~ Tb(v)/t(v). 


Figure 2: A simplified representation of the en¬ 
ergy levels of atomic hydrogen showing the hyper¬ 
fine splitting of the ground state. Ly-a and other 
higher energy level transitions may play a role in 
bringing equilibrium between the hyperfine states. 


Radio telescopes 

Though there may be considerable varieties (see, e.g., LOFAR and ORT in Figure [3]), radio 
telescopes are generally parabolic antennas to receive and detect radio frequency electro- 
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Figure 3: (A): 100 m diametre single dish Effelsberg Radio Telescope in Germany; (B): One of 
the 30 dishes of the GMRT array; (C): A few of the 25 m dishes of the Westerbork Synthesis Radio 
Telescope (WSRT; a 14 element 2.7 km East-West array) in the Netherlands. Not all radio telescopes 
are dish shaped. (D): The Ooty Radio Telescope (ORT) is a 530 m x 30 m cylindrical paraboloid 
on equatorial mount; (E): Part of a Low Frequency Array (LOFAR) station in Germany consists of 
the low band dipole antennas. See text in the “Radio Telescope” sub-section for details. 

magnetic waves. A radio telescope may consists of either a single element (“single dish” 
telescope), or multiple elements (“array” telescope). The technique of interferometry is 
routinely used in radio astronomy to achieve better sensitivity and angular resolution (or 
the “beam”) by combining signals from widely separated elements of an array telescope. 
The resolution of a single dish telescope is ~ X/D, while the same for an interferometer 
is ~ X/B for an observing wavelength of A, dish diametre I) and maximum separation 
between antennas B (known as interferometer “baseline”). A single dish telescope is sensi¬ 
tive to the total intensity / (Z, m) towards the direction of observation (Z, m) on the sky. An 
interferometer, on the other hand, directly measure the Fourier transform of / (Z, m), 



( 1 ) 


Each combination of antenna pair sample the real and the imaginary part (or amplitude and 
phase) of the visibility function V (u, v ) at an inverse angular separation (u. v) correspond¬ 
ing to projected baseline separation perpendicular to the source direction. As an example, 
Figure |4](A) and (B) shows the calibrated visibility data for a GMRT 610 MHz observation. 
At a given time, an array with N elements sample the visibility function on A r (A r — l)/2 
points on the (u, v) plane. During the course of the observation, as the apparent position 
of the source on the sky changes with time, the projected baseline separation also changes. 
This technique of sampling the (u, v ) plane is known as the aperture synthesis. The maxi- 
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Figure 4: Calibrated visibility data and image from the GMRT 610 MHz observation of the field 
containing the bright unresolved source NVSSJ171447 — 251435. (A): Visibility amplitudes and 
(B): phases as a function of baseline separation; (C): Image of a part of the field, created by inverse 
Fourier transform of the data, shows a double-lobed radio galaxy at the top right corner, and the 
point source with ~ 5 Jy flux density at the bottom left corner. The flux density scale of the image 
is saturated at ~ 125 mJy to properly show the faint extended radio galaxy. 


mum baseline for such arrays may be from few km to tens of km, and even few thousand km 
for very long baseline interferometry (VLBI). The longest baseline of the Very Long Base¬ 
line Array (VLBA) of the National Radio Astronomy Observatory (NRAO), for example, 
is ~ 8600 km (Mauna Kea, Hawaii to St. Croix, Virgin Islands). Inverse Fourier transform 
of V ( u , v ) is used to make the “radio image” of the sky. Figure |4](C) also shows the corre¬ 
sponding image created using the Astronomical Image Processing Software (AIPS). Details 
of methods for either single dish or interferometric observations, and standard data analy- 
sis procedures are b eyond the scope of this article. One may see iThompson et al.1 (1200 lb . 
Wilson et al. (2013), and references therein, for a comprehensive discussion on these topics. 


The Giant Metrewave Radio Telescope 

The GMRT is, at present, world’s largest low frequency radio telescope. It is located near 
the village Khodad in Maharashtra, approximately 80 km from Pune. The GMRT consists 
of 30 parabolic dishes, each with a diametre of 45 m. The telescope operates at ~ 150, 
235, 325, 610 and 1420 MHz. Half of the antennas are placed within central 1 km x 1 km 
area, and the rest are placed in roughly a “Y” shaped array (with largest baseline separation 
of ~ 30 km) for optimum sampling of the {u. v) Fourier domain. The unique feature of 
GMRT, apart from its high angular resolution and low frequency capabilities, is the simul¬ 
taneous sensitivity of the array to both large and small scale structures due to this hybrid 
array design. The parabolic surface of the antennas consist of wire mesh which are good 
reflector at the operating frequency range of the GMRT. All the antennas can be moved 
to point at different directions in the sky using an altitude-azimuth mounting system. The 
signals from all the antennas are transmitted to the array control building using optical fiber 
network, and combined in the correlator system to record the visibilities in a digital format. 
The angular resolution of the GMRT at 1.4 GHz is < 2 arcsecond. It is possible to easily 
achieve a spectral resolution of < 2 kHz per channel (~ 0.4 km s“ More details on 


GMRT, including all technical details are summarized by Swarup et al.l ( 199 lb . 
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Emission and absorption spectra 

A single dish radio telescope with receiver covering the 21 cm line (“L band”) can be used 
to observe radio emission from the Galactic H I. Such observations typically have coarse 
angular resolution (tens of arcminute beam size). One may also use a single dish telescope 
to observe H I in absorption towards background continuum sources. However, to do that, 
one must carefully model the emission spectra from surrounding pointing, and subtract the 
contamination from emission within the beam to estimate the absorption spectra. Thus, 
effectively one assumes that the H I emission is smooth over the scale of a few beam size. 
This assumption may not necessarily be correct. A much reliable way to get the absorption 
spectra is to carry out interferometric observations. With the typical resolution of few arc- 
second, the diffuse H I emission is resolved out, and one gets uncontaminated absorption 
spectra towards background continuum sources. Following standard data reduction steps 
(removing bad data and radio frequency interference, calibrating instrumental and atmo¬ 
spheric effects, imaging of the visibility data using inverse Fourier transform, converting 
absorption intensity to optical depth r etc.) one finally gets both the H I emission and 
absorption spectra T B (v ) and r{y), from single dish and/or interferometric observations, 
with certain spectral resolution Av. These spectra are also often expressed as a function of 
physically more meaningful Doppler velocity ( v ), instead of observed frequency (z/). 


Derived physical quantities 

For an isothermal H I cloud, both the emission and the absorption spectra can be mod¬ 
elled as Gaussian profile along the velocity axis. The central velocity of the Gaussian 
represents the line of sight component of the “bulk” velocity, and is closely related to the 
Galactic rotation. The velocity width of the line has contribution from thermal and non- 
thermal velocity dispersion. In absence of any non-thermal motion, the line width provides 
T k = (mn/k)(a/l km s -1 ) 2 K. In reality, however, due to non-thermal broadening, it only 
provides an upper limit (T k max ) of the ki netic temperature. The HI column density (i.e. the 
integrated number density along the line of sight) is given by, 


N(R I) = 1.823 x 10 18 T s / r{v)dv « 1.823 x 10 18 / T B (v)dv 


18 


( 2 ) 


where N(K I) is in cm 2 and v is in km s 1 (IKulkarni & HeilesUl988nDickev & Eockmam 
1990) . The spin temperature, T s is given by, 


T, = 


__ / T B (v)dv 

1.823 x 10 18 f r(v)dv f r(v)dv 


N(H I) 


(3) 


where t{v) is generally determined from H I absorption spectrum towards a background 
source, and N (H I) is determined from the emission spectrum for a nearby line of sight. 

In practice, if the gas is not homogeneous (i.e. having variation of density and temper¬ 
ature, along and transverse to the line of sight), or the optical depths are large (r > 1 ), 
Equation [3] can not be used to determine ‘the’ spin temperature easily. For example, mul¬ 
tiple optically thin (r < 1 ) components along the line of sight will only give a column- 
density-weighted harmonic mean temperature of the individual components. In general, 
there may not be any unique interpretation of the data, although modelling the spectra with 
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many Gaussian component s (corresponding to multip le “clouds”) is common in such situ¬ 
ations dMebold et all 1982 : Heiles & Troland . 2003al O). For extreme situations with very 
complicated profiles, where multi-Gaussian fit may be highly non-trivial and degenerate, 
one may still extract some useful information by computing T s (v) = T B (v)/r(v ) for each 
spectral channel. Note that the analysis method outlined above can be extended easily for 
imaging widespread emission or absorption against extended background sources. In those 
cases, t(u) and T B (v) will also be functions of sky coordinates (/, m). 

Figure © shows an example of 
H I absorption and emission spec¬ 
tra towards quasar PKS 1245 — 19, 
a radio bright, unresolved extragalac- 
tic source (l = 302°, b = +43°.5 
in the Galactic coordinate system). 

The emission spectrum is taken from 
the Leiden/Argentine/Bonn Survey 
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absorption spectrum is from the 
GMRT. One can clearly see multi¬ 
ple emission and absorption compo¬ 
nents at different velocity (ulsr) with 
respect to the local standard-of-rest 
(LSR) frame. The H I column den¬ 
sity, A r (JI I), computed from the in¬ 
tegral of the emission profile, is ~ 

3.7 x 10 20 cm -2 , and the average 
spin temperature is (T s ) = 1288 ± 

41 K. Both the emission and absorp¬ 
tion spectra can also be modelled as 
a set of multiple Gaussian compo¬ 
nents. There are some broad emis¬ 
sion components without correspond¬ 
ing absorption components. Such 
components are believed arise from higher temperature gas (i.e. W NM), with opti 


Figure 5: Profiles of the observed H I 21 cm emis¬ 
sion and absorption towards l = 302°, b = +43°.5. 
The emission spectrum is taken from the LAB survey, 
and the absorption spectrum is obtained with the GMRT 
(also see the section titled “Temperature of the diffuse 
neutral ISM”). 


cal depth lower than the detection limit of the observation (IRadhakrishnan et all 1 19721: 


Heiles & Trolandl 12003 aUbh . For this particular line of sight, from (T s ), and T/, lnax values 
derived from the line widths (Roy et all 2013b), one can conclude that only < 16% of H I 


is in the CNM phase (with T s < 200 K). 

Turbulence in the diffuse neutral ISM 

Indications of turbulence in the diffuse neutral ISM come mainly from observations of 
small scale fluctuations of H I emission as well as of optical depth. Recently, it has been 
clearly shown from an ongoing H I absorption survey dRov et al. . 2013a b!, also see the 
section titled “Temperature of the diffuse neutral ISM”) that, for the diffuse H I, the non- 
thermal broadening contributes significantly (i.e. comparable to the thermal broadening) to 
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the observed line width. There are direct observational evidences of AU scale (milliarcsec 
angular scale) structure s in H I opacity flu ctuations from VLBI observations (IBrogan et all 


2005: iRov et all 2012; Dutt a et all 12014). and also from variation of absorption towards 


high proper motion pulsars detected using multi-epoch observations (IFrail et all. 1994). 


Deshpande et al. (2000) reported scale-free opacity fluctuations at pc and sub-pc scales 
using interferometric observations towards extended sources. Similarly, H I emission stud¬ 
ies have revealed scale-free structures of the neutral ISM in the Mi l ky Way and other nearby 
galaxies ( Crovisier & Dickey . 1983 : Green . 1993 : Westpfahl et al. . 19991: Dutta etal. . 2009 ). 
Most of these studies use second order statistics (e.g. structure function, power spectrum) 
of int ensity and opacity fluctuation s , and c ompare those with analytical or numerical m od- 


els (Kolmogorov, 119411 : 


1964; Kraichna, 1965; Goldreich & Sridhar, 


to 


constrain the nature of the turbulence. 

Power spectra of intensity or opacity fluctuations are related directly to the underlying 
density fluctuations, but also got affected by the velocity fluctuations (iLazar i an & Pogosvan . 
2000). The power spectrum, which is the Fourier transform of the autocorrelation function, 
is given by 

P(u,v) = [ [ £(l,m)e- 2ni{ul+vm) dldm (4) 


where £(/, m) is the autocorrelation function of the quantity of interest (intensity, opacity, 
density or velocity fluctuations). The function £(/, m) is also closely related to the structure 
function S(l,m). For the intensity fluctuation, £ and S can be written as 


£(/ — /', m — m!) = (5I(l,m)8I(l',m')) 

S(l — l', m — m!) = ([1(1, m) — 1(1', m')] 2 ). 


(5) 


Here, angular brackets denote ensem¬ 
ble average over possible realization. 
If one assumes statistical isotropy, 
then P(u, v ) will be a function of only 
the magnitude U = \/u 2 + v 2 . It 
is important to note that, due to the 
scale free nature of the turbulence, the 
power spectrum is expected to have a 
power law form P(U) = AU a over 
the inertial range of scales. The struc¬ 
ture function then will also be a power 
law ( Lee & Jokipii , 1975h , with power 



law index /3 = — (a + 2). 

One advantage of interferometric 
observations is that the observed vis¬ 
ibilities are already measurements of 
Fourier transform components of the 
sky brightness. So, one can estimate 
the power spectrum from the observed 
visibilities without g oing into the complications and limitations of the imaging and decon 


Figure 6: Estimated intensity fluctuations power spec¬ 
trum for one spectral channel using visibility data from 
the GMRT observation of HI absorption towards Cas A. 
The solid line is the best-fit model with a power law 
power spectrum for the opacity fluctuations. 


volution algorithms ( Crovisier & Dickevl 1 19831: iGreenl. 1 19931: Bharadwaj & Sethi . 12001 1 
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Cas A GMRT L-band 


Based on this principle, a suitable method has been developed for radio interferometers 
like the GMRT dRov et al. . 2009 ). and is used to estimate, directly from the visibility data, 
the H I opacity fluctuations power spectrum towards the extended background emission 
from the supernova remnant Cassiopeia A (Cas A). Figure [6] shows the derived power 
spectrum for one of the velocity channels, and the best fit model overlaid on the data. 
It was established that the H I opacity fluctuations have a power law power spectrum with 
a = —2.86 ± 0.10(3cr) over the scale of ~ 500 AU to 2.5 pc ( Rov et al. . 2010n . The line of 
sight towards Cas A probes deep absorption from the local arm and the Perseus arm of the 
Milky Way. There was no significant difference in a for these two spiral arms. The results 
are in go od agreement with other studies suggesting a ~ —2 . 7 to —3 for t he d iffuse neutral 
medium ( Crovisier & Dickev . 1983 : Deshpande et al. . 2000 : Dutta et all l2009h . Although 
the power law behaviour suggests a Kolmogorov-like turbulence, the power law index is 
significantly_shallower than the Kolmogorov —11/3 spectrum for incompressible turbu¬ 
lence (IKolmogorov , 1941). This mismatch may be due to the fact that even small ionization 
fraction can couple the gas with the interstellar magnetic field. In that case, the turbulence 
wifi be magnetohydrodynamic (MHD) in nature, and the density fluctuation power spec¬ 
trum may have a different power law index (see Rov et al] 20091 and references therein for 
a detailed discussion). 

Apart from the density fluctu¬ 
ations, these data were also used 
to constrain the velocity fluctua¬ 
tions usin g “vel ocity channel anal¬ 
ysis” ( Lazarian & Pogosvanl 2000l) . 

However, the effect of the velocity 
width of the spectral channels was 
found to be insignificant on the power 
law index. This result only put a 
weak constrain on the velocity struc¬ 
ture function index /3 = 0.2 ± 0.6, 
broadly consistent with Kolmogorov 
( 1941 ) prediction of (3 — 2/3. 

It is worth mentioning here that 
this study also revealed the presence 
of MHD turbulence in the supernova 
remnant from the continuum image of 
the background source Cas A. This is 

revealed by the power spectra of the synchrotron radiation intensity fluctuations derived 
from the GMRT 1.4 GHz observation of Cas A (see Figure [7]), and also from the VLA 5 
GHz archival data for Cas A and Crab Nebula. For Crab Nebula, the power law index is 
3.23 over the whole range of U probed in these studies. For Cas A, the power spectrum 
has the same slope at large U. However, below ~ 10 kX, there is a break in the spectrum, 
and the index changes to —2.22. This break and flattening below 10 kX (25 arcsec; an an¬ 
gular scale that corresponds to the shell thickness of Cas A) ^interpreted as a 3-D to 2-D 
transition of the turbulence in the supernova shell ( Rov et all 20091) . 

Another way to get a handle on the nature of turbulence is to study the scaling of turbu¬ 
lent velocity dispersion with size of clouds. Similar studies for the molecular clouds were 
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Figure 7: Power spectrum of the synchrotron emis¬ 
sion intensity fluctuations for Cas A at 1.4 GHz esti¬ 
mated from the GMRT observation. The power law in¬ 
dex changes from —2.22 to —3.23 at U > 10 kX. 
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earlier reported by Larson (1981) and others. As mentioned before, for the dense CNM 
phase, T s is a good proxy for T k (see sub-sections titled “H I 21 cm radiation” and “Derived 
physical quantities” for details). So, from the H I emission-absorption studies, one can es¬ 
timate the non-thermal velocity dispersion v^ T oc (T/,. rnax — T k ) ~ (T fcjmax — T s ). Most of 
the time, it is not straightforward to estimate the distance of the cloud, and hence the size, 
directly. But, if the thermal pressure P = nkT k is assumed to be more or less constant, then 
the typical length scale of a cloud 


N{ H l)/n ~ N(R I )/{P/kT k ) oc N(U l)T s 


( 6 ) 


This reasoning was used to de¬ 
rive the scaling of the form i; NT ~ 
L 13 using the published millennium 
Arecibo 21-cm absorption-line survey 


results (IHeiles & TrolandL 2003a,b). 


The results are summarized in Fig¬ 
ure [8j Over almost two orders of 
magnitude length scale, the data is 
well-represented by [3 = 0.7 ± 0.1 


(Roy et all 2008). For )V(H1)T. 


< 


1Q21.5 —2 


cm ~K, the spread is higher, 
plausibly due to systematic effects of 
the spectral and spatial resolution of 
the Arecibo telescope. The result is 
consistent with the 2/3 scalin g pre¬ 



dicted by Kolmogorov (1941), who 


Figure 8: Scaling of non-thermal line width with 
length scale derived from the millennium Arecibo 21- 
cm absorption-line survey results. A power law model 
with power index 2/3 is shown as a solid line. 


derived this relation by assuming that 
the dissipation rate is scale invariant 

over the inertial range. Note that this correlation, whether or not actually related to a 
Kolmogorov-like turbulence, can be used to estimate the non-thermal velocity dispersion 
and, in turn, the true kinetic temperature, taking it as a purely phenomenological model. 

Temperature of the diffuse neutral ISM 

In the classical model of multiphase ISM, the balance of heating and cooling rate leads 
to the two stable “phases” (CNM and WNM) of the neutral ISM. Gas at the intermedi¬ 
ate kinetic temperature ( T k ~ 500 — 5000 K) will suffer from runaway heating or cool¬ 
ing to evolve into either the WNM or the CNM respectively. So, in this model, H I 
at intermediate temperature exists “only as a transient phenomenon”, and does not con¬ 
tribute significantly to the total column density dField et al. . 19691: McKee & Ostrikeil 1977 


Wolfire et all l2003h . From an observational point of view, H I emission-absorption stud¬ 
ies have shown the presence of narrow CN M components with high opacity, and with in¬ 
ferred T k in the expecte d range of < 200 K (IRadhakrishnan et all 1 19721: IHeiles & TrolandL 
2003a : iRov et al., 2006). The wide emission components with very little or no correspond¬ 
ing absorption are assumed to be originating from the WNM (see Figure [5] and the re¬ 
lated discussion in the “Derived physical quantities” sub-section). Due to their low optical 


depth, measurements of the WNM spin or kinetic temperature are rare (ICarilli et all 19 98 
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Dwar akanat h et all 120021: iHeilesL 120011: iBegum et al.l. 1201(1: Murray et all 120141) . How¬ 


ever, the thermal steady state scenario was considered to be the correct model broadly in 
agreement with the data. 

Over the last few years, there were observational indications of the presence of a sig¬ 
nificant amount of H I with temperature in the unstable range. These results seriously chal¬ 
lenged the status of the standard mode l. Unfortunate l y, only a few of these measurem ents 
are from interferometric observations ( Kanekar et al. . 2003 : Braun & Kanekaii 2005 ). and 
the rest comes from single dish observations (Millennium Arecibo 21-cm absorption-line 
survey; Heiles & Troland, 2003b). The main issue with such single dish observations is that 
the effect of the H I emission must be modelled accurately to derive the absorption spectra. 
This, and the contamination of the spectra by unwanted emission from the main beam may 
have significant systematic effects, and make the results less reliable. The other potential 
issue is that, for WNM, neither T s nor T fc max is a reliable proxy for T k ; the non-thermal 
broadening may be significant, and also T s may considerably differ from T k ( Liszt! 2001 ). 
Finally, in presence of CNM components along the line of sight, “self-absorption” of the 
background emission by the foreground components may make modelling the emission 
spectrum, if not impossible, really difficult. Even considering these possible observational 
shortcomings, the main result, that a large fraction of gas is in the “unstable” phase, has 
serious implications for the thermal steady state model. 

It is important to note that the 


Heiles & Troland! (2003b) draw the 
conclusion of a significant fraction 
of unstable gas based on T fc max de¬ 
rived from the observed emission 
line widths. It is possible to use 
a Kolmogorov-like scaling (see Fig¬ 
ure [HI) to estimate the turbulent ve¬ 
locity dispersion, and derive the true 
temperature distribution. A reanalysis 
of the Heiles & Troland ( 2003b ) data 
implementing this correction notice¬ 
ably alters the distribution ( Roy et al .. 


i 

z 



T < 500 K 


500 K < T < 5000 K 
Temperature 


T > 5000 K 


2009). Figure [9] shows the main result 


- the distribution of N (H I) at differ¬ 
ent temperature ranges before and af¬ 
ter correcting the non-thermal broad¬ 
ening. As expected, a lot of “warm” 
gases based on 7% max are indeed just 
cold gas with significant non-thermal 
broadening. There is, in fact, very lit- 


Figure 9: Distribution of H I in different phases based 
on temperatures computed without and with correcting 
the non-thermal broadening (bars marked with T),. max 
and T k respectively). Red (blue) arc H I emission com¬ 
ponents that are detected (not detected) in absorption. 
Note that the correction implies only a small amount 
of gas in the stable WNM phase; however, the unstable 
fraction remains almost unchanged. 


tie (« 10%) H I in the stable warm 

phase in this sample. However, even after this correction, a similar fraction of gas (« 30%) 
remains in the unstable phase (I Rov et al .. 2009). 

The fact that, due to various systematic effects, the single dish absorption spectra re¬ 
mains less reliable (particularly for weak and wide WNM components), prompted us to also 
start an interferometric survey with the GMRT and the WSRT. The main aim of this ongoing 


12 















































Figure 10: (A): H I 21 cm optical depth at 1 km s -1 resolution from the absoiption survey, and 
the corresponding Te from the LAB survey. Solid lines mark different T s values. Channels with 
intermediate T s may be a mix of cold and warm gas. (B): Integrated column density and average 
spin temperature for all 34 lines of sight. The high T s values imply presence of warm gas. 


survey is to acquire high sensitivity, high spectral resolution Galactic H I 21 cm absorption 
towards a large number of background sources. A large sample, and a high sensitivity 
(RMS opacity ~ 10 -3 ) are necessary to draw a statistically significant conclusion about the 
WNM properties. So far, ~ 50 lines of sight have been observed (6 — 8 hr time per source). 
Details of the survey, the analysis procedure and results based on a brighter subsample of 34 
sources, including two southern sources observed with the Australian Telescope Compact 
Array (ATCA), have been reported recently dRov et al. . 2013alibi) . One example absorption 
spectrum from this survey is shown in Figured 

From each of these lines of sight, using absorption spectra from this survey and the cor¬ 
responding emission spectra from the LAB survey, T s values are computed at ~ 1 km s ' 1 
velocity resolution for each channel with > 3 cr detection of H I absorption. As shown in 
Figure [10] (A), many channels have T s > 200 K, clearly an indication of a mix of cold and 
warm gas. The line of sight average T s , plotted in Figure [TO] (B), also have values > 200 
K for most of the cases. Based on T s at 1 km s _1 resolution, about ~ 50% of the gas (by 
column density) has T s < 200 K (i.e. channels with large cold gas fraction), and > 10% 
has T s < 1000 K (possibly with large warm gas fraction). The rest, ~ 40% of the gas, 
with 200 K < T s < 1000 K is not necessarily “unstable” gas, but can be only a mix of 
CNM and WNM. However, even without Gaussian decomposition of the profiles, based 
on the average T s values and the sensitivity limit of the survey, one can conclude that the 
detected H I absorption is not entirely from cold gas with T k < 200 K, and must also ha ve 
contribution from (stable or unstable) warm gas (Kanekar et al., 201 1 : IRov et al .. 2013al b). 
An upper limit of the CNM fraction was derived assuming that all cold gas has T s = 200 
K. The median CNM fraction was found to be / CNM ~ 0.52 for the sample. 

Finally, Gaussian decomposition of the absorption profiles were carried out, using the 
Levenberg-Marquardt ^-minimization, to estimate TT max of the individual clouds along 
the lines of sight. It was found that the number of components with T k max > 5000 K is 
very few (only seven in this range, and, considering the measurement uncertainty, possi¬ 
bly a few more out of a total 214 components). As the survey has adequate sensitivity to 
detect H I in absorption even at this temperature range, one can put a strong constraint on 
the fraction of gas in the “classical” warm phase. Combining the fact (i) that only a few 
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components have T fcimax in the range for stable WNM, (ii) that a significant fraction of gas 
can not be below the detection limit of the survey, and (iii) that the median /cnm is at most 
about 0.5, the lower limit of unstable gas fraction was found to be ~ 0.3. We conclude that, 
indeed, a significant fraction (> 30%) of the gas is in the so-called unstable phase. 


Discussions and Conclusions 


Turbu lenc e is an important ingredien t of ISM physics in various regards ('see lLlmegreen & Scalo 
20041 and lScalo & Elmegreenl 2004. for comprehensive reviews on observations and impli¬ 
cations of the ISM turbulence). The small scale structures of the ISM, originating from 
turbulence, seems to be ubiquitous not only in the Milky Way, but also in other galaxie s 


(IWestpfahl et all 1 19991: lElmegreen et all 1200 ll : iDutta et all 2009; Srianand et all 120131) . 


Fro m the galaxy disk scale to molecular cloud, i t affec t s the ISM dynamics and pr oper- 
ties ( Romeo et al. . 2010 : Shadmehri & Elmegreen . 2011 : Hoffmann & Romeo . 2012 1. The 
results reported in this article, based on H I observations, show that both the density and 
velocity fluctuation statistics are consistent with a near-Kolmogorov turbulence in the dif¬ 
fuse ISM of our Galaxy. A working model is also proposed to estimate the non-thermal 
broadening of the line width and to derive the true kinetic temperature of the gas. 

Apart from the H I 21 cm line, various other tracers at different wavelengths are also 


used to study turbu lence in different phases of the ISM (We ltv et all! 1994 iGibsonl 12007 


Boisse et all 2013 ). However, neutral hydrogen, which is present almost everywhere in 
the Galaxy, can be used to probe turbulence over a wide range of scales, using multiple 
observational techniques (from single dish emission study to VLBI small scale absorption 
measurement). For example, VLBI observations can constrain the density and the velocity 
power spectra at milliarcsec scales ( Rov et al. . 20121: IDutta et al. . [20141) . Constraining the 


turbulent behaviour of the ISM at the scale of VLBI observation can, in near future, address 
whether the same power law behaviour of the energy cascade is observed at pc and AU 
scales, and also constrain the energy dissipation scale and mechanism. It will also be useful 
to see, how the VLBI observations at AU scale compare with observations of tiny scale 
structures detected using multiwavelength observations of other tracers . 

Recent Galactic H I observations are also questioning the classical thermal steady state 
model of the multiphase ISM. There is clear evidence that > 30% of the diffuse neutral 
gas is in the unstable phase. However, any effort of understanding the properties of the 
multiphase ISM should also take into account the possible effects of the turbulence. Obser- 
vationally, one must get a handle on the turbulent velocity dispersion to accurately estimate 
the temperature, and to critically examine the validity of thermal steady state model. From 
theoretical consideration, on the other hand, turbulence may potentially play a role in cre¬ 
ating and supporting gas at the intermediate temperature range. Recent numerical studies 
suggest that turbulence can explain the presence of some amount of gas in the unstable 
warm phase (lAudit & Hennehelld. I2005L120101) . These numerical simulations, however, 
still predict a large fraction of gas to be in the stable WNM phase, somewhat in contrast 
to the observational results. More recent hydrodynamic simulations by Kim et al. (2014) 
can explain some of the observational results of Roy et al. ( 2013allbl) . It remains interesting 
to see whether more data and further numerical studies will result in a better agreement 
between the details of the observational results and the theoretical predictions. 
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